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The solid immersion lens is a powerful optical tool that allows light entering material from air or 
vacuum to focus to a spot much smaller than the free-space wavelength. Conventionally, however, 
they rely on semispherical topographies and are non-planar and bulky, which limits their integra¬ 
tion in many applications. Recently, there has been considerable interest in using planar structures, 
referred to as metasurfaces, to construct hat optical components for manipulating light in unusual 
ways. Here, we propose and demonstrate the concept of a planar immersion lens based on meta¬ 
surfaces. The resulting planar device, when placed near an interface between air and dielectric 
material, can focus electromagnetic radiation incident from air to a spot in material smaller than 
the free-space wavelength. As an experimental demonstration, we fabricate an ultrathin and hexi- 
ble microwave lens and further show that it achieves wireless energy transfer in material mimicking 
biological tissue. 
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I. INTRODUCTION 

When light is focused from air into material, refraction at the air-material interface determines the diffraction 
limit. Conventional optical lenses, placed in the far-held of the interface, control only propagating wave components 
in air. As a result, their focusing resolution in material is diffraction-limited at the free-space wavelength A because 
higher wavevector components in material cannot be accessed by far-held light [l| . These high wavevector components 
correspond to plane waves propagating at angles greater than the critical angle, which are trapped in the material by 
total internal rehection [Fig. [T]^a)]. 

Solid immersion lenses are semispherical domes of high-index material placed at or near the air-material interface 
that allow light to access these “forbidden” angles of refraction. This capability enables light to be focused to a spot 
much smaller than the free-space wavelength, with a dihraction-limited resolution set by the material wavelen^h 
^ \/n [Fig. [Hb)]. Solid immersion lenses have found extensive use in many applications, including imaging [2|, 
data storage 0, and lithography Q. They are, however, intrinsically three-dimensional and bulky [5[. Replacing 
conventional solid immersion lenses with flat counterparts would afford opportunities for integration in complex 
systems, including nanophotonic chips or, in the low frequency regime, conformal biomedical devices. 

In this paper, we propose and demonstrate a planar immersion lens based on metasurfaces. Metasurfaces are flat 
devices consisting of structured arrays of subwavelength apertures or scatterers that provide an abrupt change in 
electromagnetic properties as light propagates across the surface [6|. The properties of metasurfaces can be tuned 
by varying the parameters of the individual subwavelength elements to form a desired spatially varying response. 
This extraordinary freedom in design has been used to create devices that generate negative refraction [3, p, light 
vortices [9[, flat lensing BEil, holograms BEl, and other unusual interface phenomena in both optical [l^ 
microwave dilli regimes. We report metasurfaces based on electrically thin metallic strips with deep subwavelength 
spacing that allow radiation incident from air to refract into forbidden angles in material [Fig. [He)]. We then use 
this capability to develop a thin and planar device that reproduces the functionality of a solid immersion lens. The 
device is fabricated on a flexible substrate and demonstrated at microwave frequencies. 


II. REFRACTION AT FORBIDDEN ANGLES 

To allow for interface phenomena different from classical reflection and refraction, a metasurface can be placed at 
the air-material boundary to break translational symmetry at the interface. When the metasurface imparts a phase 
with constant gradient on incident light, propagation is governed by a generalized form of Snell’s law Q- The 
law implies that radiation incident at an angle ^inc refracts at a forbidden angle |^ref| > ^critical if the phase gradient 
is sufficiently large |V^| > ko — ko sin \ 0inc\ (see Methods). 

We implement a phase gradient by non-periodically modulating the surface with subwavelength structures of varying 
impedances. Fig. [2][a) shows the metasurface consisting of metallic strips loaded with passive lumped elements 
(resistors, capacitors, and inductors). At microwave frequencies, these elements can consist of patterned metal traces 
or commercial impedance components. Across a resonance, the phase of the current in the strips differs from that of 
the driving electric field by a value between 0 and tt. By selecting suitable passive elements, taking into account both 
the intrinsic and mutual impedances of the structures, the spatial phase profile of the transmitted wave can be shaped 
within this phase range. The use of discrete passive elements considerably simplifies the design as the metasurface can 
be reconfigured by simply changing the elements [illil. Because coupling is explicitly accounted for in the design, 
the inter-element spacing can be made deeply subwavelength. The phase range can be extended to the full 0 to 27r 
by exploiting changes in polarization (Berry phase) Q, incorporating elements with a magnetic response 0113, or 
cascading multiple layers jllj, although the limited range achieved here using a single layer is sufficient to realize a 
planar immersion lens. 

Refraction at a forbidden angle is shown in Fig. [2][b) using the metasurface to create a phase gradient of = 
7r/0.55A for radiation at 1.5 GHz. The spacing between the elements is A/20, and is necessarily subwavelength in order 
to satisfy sampling requirements. For an s-polarized plane wave incident at ^inc = 30°, the beam is entirely refracted 
to an anomalous angle of Oref = 45° that lies well beyond the critical angle ^critical = 30°. Because the metasurface 
does not rely on polarization conversion, there is no co-polarized component refracted at an angle dictated by the 
standard Snell’s law 0111. Unlike diffraction gratings, which alter the spatial amplitude profile, the metasurface 
refracts the incident wave by modulating its phase profile and thus does not result in unwanted diffractive orders. 

We find that anomalous refraction still occurs when a A/20 air gap is introduced between the metasurface and 
the interface [Fig. Efb)]. This effect is closely related to frustrated total internal reflection. In absence of material 
below the metasurface, the incident wave completely reflects off the metasurface and forms an evanescent wave at the 
surface [Fig. [2][c)]. This evanescent wave propagates along the surface in the direction of the phase gradient, which is 
a non-trivial behavior that cannot be realized with a grating 0 . When the material is placed in close proximity to 
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the metasurface, the evanescent field phase matches to a propagating wave in the material, allowing the incident beam 
to tunnel into the material with the net transport of energy across the interface. By varying the angle of incidence 
relative to the phase gradient, the angular spectrum of the transmitted beam can be made to lie almost entirely in 
the forbidden region [Fig. [2][d)]. The results are in agreement with the generalized Snell’s law. The spread around 
the predicted angle is due primarily to finite size of the aperture (see Supplemental Materials). 


III. OPTIMAL FOCUSING 

In order to design a planar immersion lens, a field source in air that focuses to a \/n spot in material must first 
be found. Although focusing across a planar interface has been previously studied, classic expressions for the optimal 
field source consider only far-held light and yield a ^ A focal spot [ 22 M 23 . We introduced a more general approach in 
ref. (21) that accounts for evanescent waves at the interface. We start by formulating an optimization problem over 
the space of current sheets js in the source plane (taken to be ^ = 0). The solution to the problem is dehned to be the 
current sheet that maximizes a metric for the degree of focus [25[. For present purposes, we assume that the material 
is dissipative, allowing small but non-zero loss. We then take the efficiency of work performed on the material as the 
focusing metric 


_ a"\E{rf)\^ 

f dr e" |E(r)p 


( 1 ) 


where rj is the focal point, a" the imaginary part of the polarizability of the object at the focal point, and e" the 
imaginary part of the material dielectric permittivity. We set a to be the polarizability of a “virtual” sphere centered 
at the focal point: the sphere has the same dielectric permittivity as the background material and can be made 
arbitrarily small (e.g., the diameter of a computational mesh unit). The electric field E can be found by propagation 
from the current sheet as described by the Green’s function G(r, r'). 

The optimization problem can now be considered in the operator formalism. Using Dirac bra-ket notation, we 
represent E and respectively as functions j^/;) and \(j)) in Hilbert space. They are related through the operator 
expression l?/^) = G\(j)) where G is the Green’s function operator. We also define a focal position operator 4> such that 
the numerator in Eq. (1) can be written as 


= a" j dr (5(r/)|E(r)|2 = a" |E(r/)|^ (2) 

Similarly, a power loss operator E can be defined to yield 

(t/>|S|t/>) = Jdr e" |E(r)|2. (3) 

Optimal focusing occurs when the choice of the source current density |0) maximizes Eq. (1). Eocusing can thus be 
posed as the optimization problem 


maximize 


{cP\G^tG\4>) 


(4) 


where S is the set of all current sheets on a plane above the z = 0 plane. The form of Eq. o is a generalized eigenvalue 
problem involving the operators A := and B := G^EG. The solution is given by the two-dimensional current 

density that satisfies A|0opt) = Amax^l^opt) where A^ax is the largest generalized eigenvalue. If B is invertible, then 
the solution |0max) can be obtained from a standard eigenvalue decomposition of the operator B~^A. Numerical 
computation can be considerably accelerated by (i) selecting the plane wave basis, which diagonalizes the Green’s 
function operator for the multilayer geometry, and (ii) exploiting degeneracies due to azimuthal symmetry about the 
focal axis. The calculation reduces to inversion of dyads at each spatial frequency, without need to explicitly form 
the full system matrices. This inverse filtering process is closely related to time-reversal [26| and can be generalized 
to transparent media by allowing the material loss to asymptotically approach zero. 

We first consider the two-dimensional geometry in Eig. [3] where the material has a refractive index n = 2. Eor 
incident s-polarized radiation, we numerically solve Eq. m to obtain a source that focuses to a line at a 4A/n 
distance. A linear metasurface is used to shape a normally incident plane wave such that the exiting field matches 
the solution. The required impedance values of the passive elements are solved by using a point-matching method 
(see Appendix C). Eig.[3][c) shows that line width of the focal spot is subwavelength 0.42A EWHM. To verify that the 
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focusing effect is due to phase (not amplitude) modulation of incident wave, the passive elements are removed such 
that the surface acts as a grated aperture. The focal spot for the grating is not subwavelength A); the intensity at 
the focal point is also decreased by a factor of 4 [Fig. ^c)]. 

The physics underlying the leasing effect is substantially different from near-held focusing devices. Unlike near-held 
plates, which focus evanescent waves at a strictly subwavelength distance (typically less than A/10), our lens’ focusing 
ability results from conventional interference between propagating waves and, as a result, the focal plane can be many 
wavelengths away. The enhanced resolution of our lens follows from the shaping of the near-held phase prohle, which 
couples the incident wave to forbidden angles on interaction with material, rather than the near-held interference 
ehects of near-held plates. Because the focusing is not subject to the intrinsic decay of the near-held, the intensity 
at the focal spot can be comparable to or higher than the incident intensity. As with solid immersion lenses, the 
focusing resolution remains subject to the dihraction limit, although with spot size set by the material rather than 
the free-space wavelength. 

Next, we consider a three-dimensional geometry where a planar source is positioned a subwavelength distance 
(A/15) above a material whose refractive index at microwave frequencies approximates biological tissue (real part 
n = 8.8). Due to symmetry about the focal axis, the polarization of the helds at the focal point can be arbitrarily 
specihed. Setting the electric field to be linearly polarized in the x direction, the solution to Eq. is found to be a 
surface wave consisting of concentric ring-like currents around the focal axis [Fig. |4l[a)]. In air, the resulting fields are 
evanescent and non-stationary, propagating in-plane towards the focal axis (see Supplemental Video). Importantly, 
the intensity profile at the source plane is significantly non-zero only within a finite circular region [Fig. SJ^b)]. The 
radius of this region defines an effective aperture size that is directly related to the loss in the material system and 
the depth of focus. At the focal plane in the material, the field converges to a spot of width A/I I, measured full-width 
half-maximum (FWHM), at a distance of about 2.3A/n (wavelength in material) from the source plane [Fig. lU^b)]. 
Although the wave originates in air, the spot size approaches Abbe’s diffraction limit A/(2n sin6>ap) in homogenous 
material, where ^ap is the half-angle the aperture subtends the focal point, due to the source’s ability to access 
forbidden wave components [Fig. [4](c)]. 


IV. EXPERIMENTAL DEMONSTRATION 

Based on the solution to the optimization problem, we designed a planar immersion lens capable of focusing an 
incident plane wave. The lens consists of concentric rings loaded with passive elements [Fig.[5][a)] forming an aperture 
that spans the width of the calculated intensity distribution at the source plane. In agreement with time-reversal 
considerations [26|, , the phase response of the lens is parabolic with concavity reversed from that expected of an 

outward propagating wave [Fig. [5][b)]. Fig. [5][c) shows that the lens focuses to a spot of size of A/8 FWHM. The 
spot size is slightly larger compared to A/11 for the optimal current sheet due to the mismatch in the phase response 
outside of the finite extent of the aperture [Fig. [5][b)]. 

In the microwave regime, an application for the immersion lens is energy transfer through biological tissue [28|. 
Since tissue exhibits a relatively large refractive indices at gigahertz frequencies (e.g., real part of 7.4 for muscle at 
1.6 GHz) [^, the efficiency of energy transfer can be substantially enhanced if the focal spot size can approach 
^ X/n. In this context, the planar nature of metasurfaces is a key advantage because it enables fabrication on flexible 
substrates. The currents in the rings can be set up by applying an electric dipole moment across the center element, 
similar to recently reported near-held plates 0, such that the lens can be excited through a compact electronic 
source rather than an incident plane wave. 

We designed the lens shown in Fig. [6l[a) consisting of modihed concentric rings fabricated on an ultrathin, hexible 
FR4 microwave substrate (100 /im thickness) [Fig. [6][b)]. The reactive elements are solved assuming an electric dipole 
excitation across the center element. Energy transfer is demonstrated by operating the lens above a simulated tissue 
volume. The setup consists of a glass tank hlled with saline solution with dimensions large enough to minimize 
rehection. We positioned the lens I cm above the air-solution interface and scanned a magnetic held probe through 
the volume. The lens couples energy into the tank with a 90% efficiency, estimated from the scattering parameters 
(see Appendix A), and produces a highly symmetric focal spot 5 cm from the lens. Fig.[6l[c) shows the magnetic held 
intensity prohle of the focal spot. The spot diameter is about 1.7 cm (that is, A/II) FWHM and is in close agreement 
with the theoretical optimal. 

Energy at the focal point can be harvested with a subwavelength coil. We attached a I-mm radius coil on a 
microelectronic implant and tested the device in a liquid volume simulating biological tissue. We set the output 
power to 500 mW (about the power radiated by cell phones), such that radiation exposure levels are well below safety 
thresholds for humans [32|. Fig. [Sfd) shows the experimental characteristics of the system. At the focal plane, the 
transferred power is measured to be about 290 /iW. This level exceeds the power consumption of many classes of 
bioelectronics (a pacemaker, by comparison, consumes 8 /iW). The transfer depth is approximately 50 times the radius 
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of the coil. At this range, the transfer efficiency of the lens (about 6 x 10“^) is significantly higher than standard 
inductive coupling systems, which are limited to efficiencies on the order of 10“^ in similar configurations [25|. The 
enhanced performance is key for wireless powering to be practical, since the output power is usually limited (typically 
from 1 to 5 W) due to safety considerations psl ]. 


V. CONCLUSION 

We have demonstrated a planar device based on metasurfaces with the functionality of a solid immersion lens. 
The enhanced focusing resolution of the device results from the ability of metasurfaces to control the near-held with 
subwavelength resolution on interaction with dielectric material. At optical frequencies, planar immersion lenses could 
be implemented with closely-spaced plasmonic antennas Q or dielectric resonators [2l| , with mutual interactions ac¬ 
counted for by tuning the properties of optical “lumped” elements [ssl, • By incorporating subwavelength structures 
that interact with the magnetic held component of incident radiation, the metasurface could also modify the optical 
impedance, allowing rehection at the interface to be eliminated As the fabrication of the metasurface is simple 

and planar in nature, the metasurface-based lens can be integrated into complex systems, such as nanophotonic chips 
or conformal biomedical devices. 


APPENDIX A: METHODS 

Metasurface refraction. For the refraction simulation in Fig. 2, the metasurface is modeled as a linear array of 
15 copper strips in FDTD software (CST Microwave Suite). The strips are 0.5 cm wide and spaced 1 cm apart on 
a FR4 substrate of thickness 1.2 mm. Rehective metal sheets were placed around the aperture in order to remove 
edge dihraction. The structure is excited by a plane wave at varying angles of incidence at 1.5 GHz with perfect 
electric conductor boundaries on the y = ±1 cm planes. The predicted angles of refraction are calculated using the 
generalized Snell’s law 


sin ^ref — 


- sin6>inc 
n 



(Al) 


where n is the index of refraction, O^ei the angle of refraction, 6>inc the angle of incidence, and V4> the phase gradi¬ 
ent 0. The passive components are used to create a phase gradient ofV4> = 7r/0.55A; their values are given in the 
Supplemental Materials. 

Metasurface focusing. The focusing metasurface in Fig. [5] is excited by a plane wave at 1.6 GHz y = ±5 cm at 
normal incidence. The concentric rings are placed in a circular aperture of diameter 8.4 cm. The largest ring has a 
diameter of 7.5 cm. To calculate the impedance values of the passive elements, a multiport simulation of the structure 
is performed. The required passive elements are then obtained by solving a set of balance equations in the impedance 
matrix (see Appendix C). The values of the passive elements are given in the Supplemental Materials. 

Experimental setup. The prototype of the planar immersion lens was fabricated on an ultrathin (100 jam) FR4 
dielectric (e^ = 4.4) substrate using commercially available surface-mount capacitors and inductors. The components 
are mounted on the substrate using standard soldering techniques. The lens is excited by a co-axial cable attached 
by a SMA connector. The connector is soldered across the center element in the lens [Fig. 5(a)], with outer and inner 
conductors separated by the gap across the center circular element. The prototype dimensions are 8 cmx8 cm. For 
the experimental measurement, a liquid volume consisting of 0.5% saline solution was used as the focusing medium. 
The complex permittivity of the liquid was measured to be = 77.2 + 15.4i using a dielectric probe. An excitation 
signal (500 mW, 1.6 GHz) was sent to the lens from the signal generator (Agilent E4436B) through a power amplifier 
with co-axial cables. With the lens surface positioned 1 cm from the simulated tissue volume, the reflection was 
measured to be Su = —9.8 dB using an Agilent E5072A VNA. Since the back-radiation (power radiated into air) 
is measured to be < 1%, the efficiency of power coupled into the solution can be calculated as I — Su ^ 90%. The 
spatial distribution field in the material is measured by scanning a wideband magnetic field probe (10 MHz to 18 GHz; 
AEMP002, AET, Inc.) through the volume using a robotic positioner (4EM500, NEC, Inc.). 


APPENDIX B: DEGENERACY AND THE MATCHED FILTER 


The operators A := and B := are symmetric under rotation about the 2 ; axis and therefore commute 

with the transverse rotation operator Rg. The generalized eigenvalue problem in Eq. (|1]) is degenerate, since if 10) 
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is a generalized eigenvector, then any rotation Re\(t)) is also a generalized eigenvector. At the focal point, Rq rotates 
the transverse polarization of the field vector. The right-handed (RHP) and left-handed (LHP) circular polarization 
states are the eigenvectors of Re and span the space of solutions. Because any transverse polarization can be generated 
by the superposition of RHP and LHP states, the x and y components of the polarization at the focal point can be 
arbitrarily fixed without loss of optimality. In this case, the generalized eigenvalue problem reduces to a matched 
filtering problem. To see this, let h be a unit vector in the direction of the polarizability vector ah. The polarization 
current at the focal point is then 


Jp = 


dt 


aSirf) ( 


fdF, 

dt 


(Bl) 


Adopting phasor notation with time dependency exp(—icjt) 


, the rate of work performed at the focal point is 




-wa"|E(r/) -fil' 


which yields the following expression for efficiency 


J dr e"(r) |Ep 


(B2) 


(B3) 


Let = f dr a"S(r — rj,r' — rj)hh be the new focal point selection operator, where hh is a dyad. Because of the 
polarization vector h, 4>' can now be written as an outer product |x)(x|. In operator notation, Eq. (jB2p becomes 




(B4) 


which can be maximized by the matched filter. It can be shown from the Cauchy-Schwarz inequality [25| that y is 
maximized when 


|0opt) = (G^EG)-'Glx) (B5) 

For a multilayer material system, G is diagonal for the angular spectrum basis. Evaluation of Eq. (jB5[) requires only 
the inversion of individual dyads at each (kx^ky) component, without need to form the full system matrices. 


APPENDIX C: SELECTION OF IMPEDANCE ELEMENTS 

We consider a multiport structure {N ports) composed of primitive elements designed to induce a desired field 
pattern at the air-material interface. The structure of the elements are selected according to the dominant mode 
of interaction (e.g. rings for magnetic response or strips for electric response) and are spaced with an appropriate 
sampling period to discretize the output field. Let Fq be a Cartesian component of the desired field, 2 : the direction 
of focus, and z = zq the interface plane. For simplicity, consider the two-dimensional case, with the geometry being 
translationally invariant along the y axis. We begin by applying a point-matching method, similar to the design 
procedure for near-held plates [l^, to solve for the desired amplitudes and phases a G at each port. To do so, a 
set of helds corresponding to the held pattern when the nth port is excited, is hrst obtained with a held solver. 
We then construct vectors /o and /i, ..., /at by sampling the helds along the interface. From these vectors, we dehne 
a transfer matrix T = [/i, / 2 , • • •, /at]- For a set of weights a at the ports, the total held at the interface can be found 
as f = Ta. The optimal weights are found by solving a least squares problem minimizing the error between f and fo, 
which has solution Uopt = {T^T)-^T^fo. 

The goal is now to hnd a set of passive impedances {^n,n • > 0, n = 1,..., N} to terminate all or a subset 

of the ports such that port amplitudes respond with Uopt under a specihc excitation (plane wave or single port). To 
determine the values of these elements, we hrst characterize the coupling between ports with the mutual impedance 
matrix Z G We then hnd the vector of voltages vs at each port under excitation. For a plane wave excitation, 

these are the voltages measured at each port under open-circuit conditions, while for an excitation at Port 1, is the 
vector Vs = ['T’o, 0,..., 0]^, where vq is the applied excitation voltage. These voltages satisfy the balance equations 


Vs Fvi = Zic 
VI = -Znic 


(Cl) 

(C2) 
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where iq is vector of currents and Zr is a diagonal matrix with diagonal entries To find the impedances Zr^^^ 

we rearrange the above to obtain 


vi = -{I + ZrZ-^) ^ZrZ-^vs (C3) 

where I is the identity matrix. Defining a := vi/vq and eo '•= vs/vq and rearranging, we find 

ZR{Z-\eQ + a)) = -a (C4) 

To force the port amplitudes to be Uopt, the passive elements are thus required to have the impedance values 


^R,n — <^opt,n/Cn 


where Cn are the entries of the vector C '•= Z ^ (ep + Uopt)- 


(C5) 
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FIG. 1. Refraction at a planar interface between air and high-index material, (a) Ordinary refraction. Light cannot be refracted 
beyond the critical angle ^critical since the corresponding beam is trapped by total internal reflection (dashed red line), (b) 
Refraction with an immersion lens. The lens enables light to both enter (black line) and escape (dashed red line) material 
at a forbidden angle. Light can thus be focused to ~ A/n, where n is the material refractive index, (c) Refraction with a 
metasurface. The metasurface changes the amplitude and phase of light, enabling it to refract at a forbidden angle. 
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FIG. 2. Refraction at forbidden angles with a metasurface, (a) Metasurface consisting of metallic strips reactively loaded 
with impedances varying in the x direction. The strips, separated by a distance of s = A/20, are periodic in the y direction 
and impart a phase gradient of 7r/0.55A on incident s-polarized radiation, (b) Magnetic field amplitude |Re(iLx)| below the 
metasurface in the air gap and material (n = 2) for a plane wave at ^inc = 30°. (c) Magnetic field in air when the material 
volume is removed, (d) Gontour plot of the angular spectrum of the refracted beam as a function of the angle of incidence. The 
solid white line is the angle of refraction predicted by the generalized Snell’s law and the solid purple line ordinary refraction. 
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FIG. 3. Focusing across an interface with a metasurface, (a) Metasurface lens above a material with refractive index n = 2. The 
strips are assumed to be periodic in the y direction, (b) Magnetic field intensity in material for incident s-polarized radiation. 
The focal plane (dashed line) is 4A/n from the lens, (c) Intensity profile at the focal plane for a plane wave at normal incidence. 
The spot width is 0.42A FWHM for the lens (dark shading) and A for the grating (light shading). The grating consists of metal 
strips without passive elements. 



FIG. 4. Optimal focusing across an interface, (a) Focusing in material with an optimized current sheet (contour plot of |Re(i/x)| 
shown). The material is 0.5% saline water (n = 8.8 + iO.9 at 1.5 GHz). The focal plane (dashed white line) is located 2.3A/n 
(wavelength in material) from the current sheet, (b) Magnetic field intensity profile and (c) angular spectrum at the (1) source, 
(2) interface, and (3) focal planes. The shaded region in (b) denotes the width of the focal spot A/ll (FWHM). Shaded regions 
in (c) show the propagating wavevectors 1 < |/cx/A^o| < 8.8 supported by material. 
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FIG. 5. Planar immersion lens with concentric rings, (a) Metasurface lens consisting of concentric rings for focusing an incident 
plane wave. The rings are loaded with passive elements. The lens is positioned A/15 above saline water, (b) Phase profile 
of the magnetic field {Hx) exiting the lens along the x axis. The phase response of the lens matches the profile dervied the 
optimization procedure within the aperture, (c) Contour plot of the magnetic field intensity. The focal plane is located 2.3A/n 
(wavelength in material) away from the lens. Simulations are performed at 1.6 GHz. 
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FIG. 6. Focusing in a dielectric liquid using a coaxial-excited lens, (a) Focusing in a 0.5% saline solution with a metasurface 
lens. The lens is excited across the center gap (white circles) at 1.6 GHz by a coaxial cable. Separation between the lens and 
interface, 1 cm; focal plane and lens, 5 cm. (b) Photograph of the lens, (c) Magnetic field intensity profile at the focal spot. 
The theory curve is calculated using the refractive index n = 8.8 + zO.9. (d) Energy transfer to a 2-mm device as a function of 
distance. Inset images show the brightness of an LED on the device. The output power is 500 mW. Scale bar, 1 cm. 



















